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Abstract 

The precise determination of the masses of the CP-odd and -even heavy Higgs 
bosons is an important part of the study of Supersymmetry and its relation with cos- 
mology through dark matter. This note presents a determination of the A mass with 
the e + e~ — > H°A° — > bbbb process for a dark matter motivated cMSSM scenario 
with Ma = 1141 GeV at CLIC. The analysis is performed with full simulation and 
reconstruction at -^=3 TeV accounting for beamstrahlung effects. SM and SUSY 
backgrounds are considered and the effect of the overlay of yy — > hadrons events 
on the signal is studied for various assumptions for the detector time-stamping ca- 
pabilities. The di-jet mass resolution is improved by applying a kinematic fit. The 
A mass can be determined with a statistical accuracy of ~ 3-5 GeV for 3 ab~' 
of statistics and to 20 bunch crossings of yy background integrated in one event, 
respectively. 



1 Introduction 



The CLIC linear collider is expected to probe physics at the multi-TeV scale, complementing 
and extending the reach of the LHC on New Physics beyond the Standard Model. If this New 
Physics is manifested by a rich spectroscopic of new particles, as several models predict, their 
study will be a key part of the CLIC physics program. In this note we study the pair production 
of heavy supersymmetric Higgs bosons, e + e~ — > H°A° — > bbbb for the parameters of bench- 
mark point K', defined in the constrained Minimal Supersymmetric Standard Model (cMSSM) 
extension in ref. 0]]. This benchmark is characterised by very heavy supersymmetric particles, 
with kinematic thresholds for pair production between 2 and 3 TeV. The H° and A bosons have 
masses of 1141 and 1137 TeV, respectively, with a natural width T = 28 GeV. The accurate de- 
termination of these masses is essential for outlining the SUSY profile and understanding the 
connection of SUSY with cosmology, through dark matter. The A mass plays a special role in 
this respect. In fact, the relation between its mass, Ma, and that of the lightest neutralino, M a, 
determines the cross section of the XiXi ~ y ^° annihilation process in the early universe, which 
could have been responsible for determining the amount of relic dark matter we observe in the 
universe today Q. Pair production of heavy neutral Higgs bosons at the linear collider has al- 
ready been studied for several scenarios with lighter particles (3JIH. This analysis is interesting 
for assessing the detector requirements in terms of ^-tagging efficiency, since the signal has four 
b jets and small production cross section, jet clustering, di-jet invariant mass reconstruction and 
detector time stamping capabilities in presence of two-photon machine-induced background. 



For this study, signal events are generated with ISASUGRA 7.69 and Pythia 6.125 [|6j. 
The observed signal cross section at 3 TeV is 0.3 fb, accounting for the CLIC luminosity spec- 
trum. The main SM background processes are generated with Pythia 6.125. In addition, 
the irreducible inclusive bbbb SM background is generated with CompHep Q at the parton 
level and subsequently hadronised with Pythia. The number of events generated for the var- 
ious processes considered here and their cross sections are given in Table [TJ Beamstrahlung 
effects are included using CALYPSO for the CLIC 2008 accelerator parameters (5J. We assume 
the beams to be unpolarised. Events are processed through full detector simulation using the 
GEANT-4-based MOKKA Q program and reconstructed with MARLlN-based (10l processors 
assuming a version of the ILD detector ifTTI . modified for physics at CLIC. A total of ~100k 
SM and SUSY events have been generated, fully simulated, reconstructed and analysed (see 
Table [[])■ Particle tracks are reconstructed by the combination of a Time Projection Chamber 
and a pixellated Si Vertex Tracker. The momentum resolution is 8p/p 2 = 2 x 10~ 5 GeV -1 and 
the impact parameter resolution Or_<i> = (5 © M v x ) Aim. The parton energy is reconstructed 



using the PandoraPFA particle flow algorithm (121 . Performances are discussed below. For 
the 77 — > hadrons background simulation two photon events are generated with the Guinea Pig 
program JT3l using the 2008 nominal CLIC beam parameters at t/s = 3 TeV where the hadronic 
background cross section is modelled following Schuler and Sjostrand |[14l . The energies of 
two colliding photons are stored with their relevant probability and passed to PYTHIA for the 
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Table 1: Summary of simulated samples 



Process 


o 


Generator 


Nb of Events 




(So) 




Generated 


H U A U 


0.3 


ISASUGRA 7.69+ 


3000 






PYTHIA 6.125 




Inclusive 


32.5 


ISASUGRA 7.69+ 


20000 


SUSY 




PYTHIA 6.125 




H + H 


0.6 


ISASUGRA 7.69+ 


3000 






PYTHIA 6.125 




w + w- 


464.9 


PYTHIA 6.125 


30000 


z°z° 


26.9 


PYTHIA 6.125 


15000 


it 


19.9 


PYTHIA 6.125 


15000 


bbbb 


0.4 


CompHEP+ 


5000 






PYTHIA 6.125 




w + w-z° 


32.8 


PYTHIA 6.125 


7500 


z°z°z° 


0.4 


PYTHIA 6.125 


2500 



generation of the hadronic events. On average, there are 3.3 yy — > hadrons events per bunch 
crossing (BX) with Myy > 3 GeV. These are passed through the same MOKKA full detector 
simulation and overlayed to the particles originating from the primary e + e~ interaction at the 
reconstruction stage, using the event overlay feature in the LCIO persistency package ifTBTl . For 
this analysis, the yy background is overlayed only to the signal H°A° events, in order to study its 
effect on the signal reconstruction and to the peaking H + H background events. These samples 
are reconstructed with different amounts of overlayed background, corresponding to different 
time-stamping performances of the detector. The remaining backgrounds have a flat di-jet in- 
variant mass distribution in the signal region and the discriminating variables adopted are not 
significantly affected by the two-photon background. The comparison of these results allows us 
to assess some requirements on the detector timing capabilities. 

3 Event Selection and Reconstruction 

At the chosen SUSY parameters, where the e + e~ — > H°A° process is close to its kinematic 
threshold, the production cross section is small. The natural signal-to-background ratios to SM 
processes such as and tt are. ~ 1/1500 and 1/60, respectively. Therefore, a set of selection 

criteria with rejection factors of <^(10 3 — 10 4 ) must be applied. There are two main sets of cuts 
which effectively reduce the SUSY and SM backgrounds. The first is based on the event shape 
and jet multiplicity, since two-fermion and gauge boson pair production events at 3 TeV are 
highly boosted, while the heavy Higgs bosons of the signal are produced with small kinetic 
energy, giving 4-jet, spherical events. The second is /^-tagging since all four jets in the signal 
sample contain b hadrons. The bbbb irreducible SM background has a cross section of only 
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0.4 fb, which makes fr-tagging the single most effective selection criterium to isolate the signal. 

In order to reject particles which are poorly reconstructed or likely to originate from under- 
lying 77 events, a set of minimal quality cuts are applied. Only particles with p, > 1 GeV are 
considered, charged particles are also required to have at least 12 hits in the tracking detectors 
and 8p/p < 1. The event selection proceeds as follows. First multi-jet hadronic events with 
little or no observed missing energy are selected. We require events to have at least 50 charged 
particles, total reconstructed energy exceeding 2.7 TeV, event thrust below 0.95, sphericity larger 
than 0.04, transverse energy exceeding 1.3 TeV and 3 < A^ ets < 6, where N- ]et& is the natural num- 
ber of jets reconstructed using the Durham clustering algorithm |[T6l with y cut = 0.0025. These 
cuts remove all the SUSY events with missing energy and the e + e~ — > ff events. 

For events fulfilling these criteria, we perform the final jet reconstruction using two different 
approaches and we study the sensitivity of the reconstruction techniques on the two-photon 
background. The first approach is to apply an exclusive reconstruction by forcing the event into 
four jets, as traditionally done in e + e~ event reconstruction. Here, we use the Durham clustering 
algorithm and we shall refer to this as the "4-jet Durham" method. The second approach uses 
a more inclusive technique. Since yy events give rise to soft jets, it appears advantageous to 
cluster jets using a method which behaves like a cone algorithm and select only the four most 
energetic jets for performing the di-jet reconstruction. We use the anti-kt clustering algorithm 
in cylindrical coordinates ifTTl . implemented in the Fast Jet package, which we have ported as 
a processor in the Marl in analysis framework. The choice of cylindrical coordinates is optimal 
since the yy — > hadrons events are forward boosted, similarly to the underlying events in pp 
collisions for which the anti-kt clustering has been optimised. For each event, we choose the 
minimum R value at which the event has exactly four jets with energy in excess of 150 GeV. 
These jets are then used for the reconstruction, while the lower energy jets, which are mostly 
due to the yy events when the background is overlayed, are disregarded. We shall refer to this 
as the "semi-inclusive anti-kt" method. The effect of ^-tagging selection is included for the 
reconstructed jets, requiring four jets to be ^-tagged. For this analysis we use a parametrised 
response and we assume that the efficiency for jets containing decays products of b quarks, Sb, is 
0.85 with a mistag probability of 0.05 and 0.01 for c and light quarks, respectively. We study the 
dependence of the result on Eb in Section 4. Then, the di-jet invariant mass is computed. Since 
there are three possible permutations for pairing the four selected jets and the H° and A bosons 
are expected to be almost degenerate in mass, we take the combination minimising the difference 
AM of the two di-jet invariant masses and require \AM\ < 190 GeV. Since the signal events are 
predominantly produced in the central region while the yy and some of the SM background 
events are forward peaked, we only accept events for which the jet with the minimum polar 
angle, 6, has | cos 0| < 0.92. Finally, events with one or more jets having an invariant mass 
within 12 GeV of the top mass are discarded, to reject it and H + H — > tbtb backgrounds. The 
number of selected events for the signal and the main background processes at each stage of 
the event selection procedure is shown in Figure [T] We study the parton energy resolution for 
jets from signal events fulfilling our selection criteria. We compare the jet energy from particle 
flow to the energy of the b hadron which is closest in energy-momentum space. We study the 

Eh—E 

distribution of the quantity — g-^- and we determine the r.m.s. of the distribution, truncated 
to include the 90% of the entries closest to the peak, RMS90. We observe a relative parton 
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Figure 1 : Number of selected di-jets for the signal and the main background processes at each 
stage of the event selection procedure 



Table 2: RMS90 relative energy resolution for b-jets in signal events reconstructed with the 4- 
jet Durham clustering and the semi-inclusive anti-kt algorithm for different amounts of 
overlayed 77 background. 



Nb of BX 


RMS 


9o/£jet 


of overlayed 77 


4-jet 


semi-inclusive 




Durham 


anti-kt 





0.113±0.002 


0.112±0.002 


20 


0.149±0.003 


0.129±0.002 


40 


0.170±0.003 


0.133±0.002 


60 


0.183±0.003 


0.140±0.002 



energy resolution RMS9o/£j e t = 0.113±0.002 and 0.112±0.002 for jets from selected signal 
H A events in absence of overlayed background using the 4-jet Durham and the semi-inclusive 
anti-kt algorithms, respectively. This resolution is due to both the performance of the particle 
flow for jet energies as high as 1.3 TeV and also the effect of neutrinos produced in semilep tonic 
heavy quark decays and escaping the detectorj^] The evolution of the RMS9o/£j e t resolution 
on the Z?-jets with the amount of overlayed 77 events for the two algorithms is summarised in 
Table El 

The Higgs mass resolution is studied by fitting the di-jet invariant mass distribution with the 



Repeating the analysis for H A events where both bosons are forced to decay into light quark and into b quarks 
without semileptonic decays with escaping neutrinos we obtain a relative parton energy resolution RMS90 = 
0.088±0.002 and RMS 90 = 0.093±0.004, respectively. 
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Figure 2: Di-jet invariant mass distribution for signal events before (open circles) and after (filled 
circles) kinematic fit using the semi-inclusive anti-kt jet clustering. 

sum of two Breit-Wigner functions, describing the A and H° signals, folded with a Gaussian 
resolution term. In this fit we set three parameters free: the A mass, Ma, the Gaussian resolution 
term, Om, and the number of signal events, N s - lgna i, while keeping the H° to A mass splitting and 
the natural widths of the bosons fixed to their values in the SUSY model. We obtain a di-jet 
invariant mass resolution Gm = (68.5±6.2) GeV with the 4-jet clustering and (76.3±6.3) GeV 
using the semi-inclusive anti-kt. 

In order to improve the di-jet mass resolution, we apply a constrained kinematic fit. We 
use the port of the PUFITC kinematic fit algorithm HH to the Marlin framework. PUFITC was 
originally developed for reconstruction in DELPHI at LEP-2 and it has been successfully 

applied for the reconstruction of simulated linear collider events at lower energies (4J. The 
kinematic fit adjusts the momenta of the four jets as pp = apM + bps + cpc, where pm is the 
jet momentum from particle flow, ps and pc are unit vectors transverse to pm and to each other 
and a, b and c free parameters. The rescaled jet momenta minimise the quantity £,■( + 

-ij- + %), where a is the expected energy loss parameter, a a the energy spread and Ob, o c 
the transverse energy spreads as in the notation used above. For this analysis, we impose the 
constraints p x = p y = and E ± \p z \ = y/s, where the last conditions accounts for beamstrahlung 
photons radiated along the beam axis. We take the values of the energy spread from simulation. 
We only accept events with a kinematic fit % 2 < 10. After the kinematic fit, the relative jet energy 
resolution RMS 90 /£j et improves to 0.094±0.002 and to 0. 103±0.002 without and with 20 BX of 
77 background, respectively. The di-jet invariant mass resolution improves by more than a factor 
of two to C M = (27.7± 4.8) GeV and (29.8±4.7) GeV using the 4-jet and the semi-inclusive anti- 
kt method, respectively (see Figure |2]). The use of a kinematic fit also mitigates the effect of 
the overlayed 77 events on the di-jet mass resolution. Since we do impose the nominal centre- 
of-mass energy, allowing for beamstrahlung, jet energies are rescaled in the fit. The Gaussian 
resolutions of the di-jet invariant masses obtained for the raw jet energy and momentum from 
the particle flow and that for the rescaled values after the kinematic fit are given in Table [3] for 
the two jet algorithms and different amounts of overlayed backgrounds. Up to about 20-30 BXs 
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Table 3: Higgs mass experimental width after kinematic fitting for different amounts of over- 
layed 77 background using the 4-jet Durham clustering and the semi-inclusive anti-kt 
algorithm. 



Nb of BX 


Kinematic Fit 


of overlayed 77 


o M jj (GeV) 




4-jet semi-incl. 





27.7± 4.8 29.8±4.7 


5 


32.0± 5.0 30.1±5.0 


20 


54.0± 8.3 34.5±6.7 


40 


72.2± 7.4 45.4±5.6 


60 


78.6±10.9 52.5±8.2 



the kinematic fit ensures a good di-jet invariant mass resolution. Above 40 BXs of background, 
corresponding to about 1.75 TeV of total energy from 77—)- hadrons events overlayed to the 
e + e~ — > H°A° event, the di-jet mass resolution quickly degrades for the 4-jet cluster while it 
remains reasonably good for the semi-inclusive anti-kt algorithm. 

4 Results 

We consider the results of the analysis of a data set corresponding to an integrated luminosity of 
3 ab _ 1 taken at %fs = 3 TeV, corresponding to the production of ~ 800 e + e~ — > H A events. The 
invariant mass distribution of the two di-jets in events selected by the criteria discussed above 
is shown in Figure [3] for various amounts of 77 background overlayed. The SM background is 
small and flat in the signal region. The only peaking background is represented by H + H — > tbtb 
events fulfilling the selection criteria. These poorly reconstructed charged Higgs decays passing 
the anti-top and jet selection criteria give di-jet masses centred at approximately the same mass 
as the H° and A Q bosons, since M H = 1148 GeV in this model, but with a broader distribution. 
We repeat the fit to the di-jet mass distribution adding a smooth background, parametrised by 
a second order polynomial and a peaking background term, parametrised by a Breit-Wigner 
distribution. We obtain the parameters of these distributions by fitting the distributions for the 
corresponding samples. The fitted number of signal, iV s ignab and background, A^kg> events and 
the Ma mass value after the kinematic fit for the various 77 overlay conditions are summarised in 
Table |4] We observe a significant advantage in using the semi-inclusive anti-kt method over the 
more traditional 4-jet clustering. While performances are comparable in absence of background, 
the statistical accuracy of the mass measurement is more stable with the amount of overlayed 
background when using the semi-inclusive anti-kt method. The A mass is measured with a 
statistical accuracy of ~ 3 GeV in absence of 77 background, which increases to ~ 5 GeV 
for 20 BXs of overlayed background. A ~ 3 a systematic upward shift of the fitted mass 
value for 20 or more bunch crossings of background, which appears for the 4-jet clustering 
method, is removed by using the semi-inclusive anti-kt clustering and the fitted mass remains 
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Figure 3: Di-jet invariant mass distribution for (a) no, (b) 20, (c) 40 and (d) 60 BXs of over- 
layed 77 background for 3 ab _1 of integrated luminosity at y/s=3 TeV. The HA signal, 
the HA events with incorrect jet pairing and the various background components are 
highlighted. 



Table 4: Number of selected background and signal di-jets and value of the A mass after kine- 
matic fitting for different amounts of overlayed 77 background for 3 ab~ 1 of integrated 
luminosity at -^=3 TeV. 



Nb. of BXs 


Nbkg 


^V signal 


M A (GeV) 


of overlayed 77 






4-jet 


semi-incl. 





76 


222±19 


1137.4±3.3 


1136.7±3.4 


5 


77 


224±20 


1144.4±3.8 


1135.9±3.7 


20 


102 


208±20 


1160.7±6.9 


1139.9±5.4 


40 


96 


183±20 


1167.2±8.2 


1134.1±7.2 


60 


97 


167±21 


1170.6±9.7 


1132.9±8.6 
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within less than 1 a from its simulated value. Finally, we repeat the analysis by changing the b- 




Figure 4: Statistical uncertainty on Ma as a function of the ^-tagging efficiency per jet, e&, for 
20 BX of 77 background overlayed. 

tagging efficiency per jet, £},, in the range 0.75 < £b < 0.95 and study the change in the statistical 
accuracy on Ma as a function e& (see Figure [4}. 

5 Conclusions 

The process e + e~ — > H°A° — > bbbb has been studied for the parameters of a dark matter moti- 
vated, high mass SUSY scenario with Ma = 1141 GeV. The analysis is based on fully simulated 
and reconstructed events and includes both physics backgrounds and machine-induced 77 events. 
CLIC at 3 TeV can study the production of pairs of heavy Higgs particles with high accuracy up 
to masses close to the production kinematic threshold. Event shape discriminating variables and 
jet flavour tagging efficiently reduce the SM and SUSY backgrounds to achieve a good signal- 
to-background ratio in the signal region at high mass. A semi-inclusive jet clustering using a 
cone-like algorithm in cylindrical coordinates helps in reducing the impact of overlayed events 
and ensure a stable, unbiased mass determination up to largest number of overlayed bunch cross- 
ings considered in this study. Due to the favourable kinematics of the signal events, the di-jet in- 
variant mass resolution can be significantly improved by performing a constrained kinematic fit, 
accounting for beamstrahlung. This fit further mitigates the impact of overlayed 77 background 
up to ~ 20-30 BXs, corresponding to detector time stamping accuracies of 10-15 ns. Under 
these conditions, the A mass can be determined with a statistical accuracy of ±3-5 GeV for an 
integrated luminosity of 3 ab _1 . This accuracy matches that required for obtaining predictions 
on the relic dark matter density in the universe from collider data with a precision comparable 
to that presently achieved in the study of the cosmic microwave background with the WMAP 
satellite data IPT91 . 

The results reported in this note can be further improved. The particle flow algorithm still 
needs to be optimised for jets of 6(\ TeV). Explicit ^-tagging has yet to be performed and can 
profit from the considerable decay length of the energetic b hadrons. The jet clustering algorithm 
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can be further optimised for the reconstruction of events in presence of a diffuse background, as 
in the case of a large number of overlayed 77 events. 
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